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The study of double vector charmonium meson production at B-factories 

within light cone formalism. 
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In this paper the processes e + e~ — > J/^J/fy, J/^ip', i/j'ip' are considered in the framework of light 
cone formalism. An important distinction of this approach in comparison to the approaches used in 
other papers is that relativistic and leading logarithmic radiative corrections to the cross section can 
be resummed within light cone formalism. In this paper the effect of this resummation is studied. 
It is shown that this effect is important especially for the production of higher charmonium mesons. 
The predicted cross sections are in agreement with the upper bounds set by Belle collaboration. 
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^ ; I. INTRODUCTION 

The measurement of the cross section of the process e + e~ — > J /^rj c at Belle collaboration [l[ revealed large 
discrepancy between the experiment and the leadingorder NRQCD prediction. Latter measurements of the processes 
of double charmonium production at B-factories [llbl shown that there is disagreement between theory 0, IE E| and 
experiment in other processes. Only in few years it was realized that the contradiction between NRQCD prediction 
and experimental result for the process e + e~ — > J[^rj c can be resolved if one takes into account radiative corrections 
i— i ! and relativistic corrections simultaneously [E Q • 

• In addition to NRQCD [nj, hard exclusive processes can studied within light cone formalism (LCF) [TTI . [l2| . Within 
Oh 1 LCF the cross section is built as an expansion over inverse powers of characteristic energy of the process. There are 
two very important advantages of LCF in comparison to NRQCD. The first one is connected with the following fact: 
LCF can be applied to study production of any meson. For instance, it is possible to study production light mesons, 
such as 7T mesons, or production heavy mesons, such as charmonium mesons. From NRQCD perspective, this implies 
that LCF resums infinite series of the relativistic corrections to amplitude, which can be very important. The second 
, advantage is that LCF resums very important part of QCD radiative corrections - the leading logarithmic radiative 
■ corrections to amplitude ~ a s log(Q). This is very important advantage since at high energies the leading logarithmic 
I/"") \ corrections can be even more important than the relativistic ones. 

The first attempts to study double charmonium production at B-factories in the framework of LCF were done 
in papers [l3l . IlH . Il5| | . The main problem of these papers is connected with rather poor knowledge of charmonium 
. . distribution amplitudes (DA). It should be noted that within LCF the calculation of hard exclusive charmonium 
production cannot be considered reliable if one has poor knowledge of DAs. Fortunately, lately charmonium DAs 
became the object of intensive study [H, HE EE EE EE El El . In papers [IE EE EH the models of DAs for 15 and 
2S states charmonium mesons were proposed. If one uses these DAs to calculate the cross sections of the processes 
e + e~ — * J/^r] c , J/^r)' c ,ip'r) c ,ip'r)' c , the agreement with the experiments can be achieved [HI]. In present paper these 
models of DAs will be used. 

A lesson that can be learnt from the study of the process e + e" — * J/^ijc within NRQCD and LCF is that the 
leading order NRQCD predictions for hard exclusive processes cannot be considered as reliable before the relativistic 
and radiative corrections are not taken into the account. This paper is devoted to the study of the hard exclusive 
processes e + e~ — > J/\[ r J/\E r , J/tyip', ip'ip* in the framework of LCF. An imp ortant distinction of this paper in 
comparison to the papers where these processes were studied earlier [24|, EE EE EH EE E3| is that within LCF the 
relativistic and leading logarithmic radiative corrections to the cross section can be resummed. Thus one can hope 
that the predictions obtained in this way are more reliable. 

This paper is organized as follows. Next section is devoted to the calculation of the cross sections of the processes 
under study at the leading order approximation of LCF. In the third section 1/s corrections to the leading order result 
will be considered. In the last section the result of the calculation will be presented and discussed. 
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FIG. 1: The diagrams that contribute to the process e + e — * Vi(pi)V2(f>2) at the leading order approximation in strong 
coupling constant. 



II. THE LEADING ORDER CONTRIBUTION. 



The diagrams that contribute to the process e + e~ — > V^,(pi)Va(p2) &t the leading order approximation in a s are 
shown in Fig. [TJ The diagrams shown in Fig. [TJi, b can be divided into two parts. The first part is the annihilation 
e + e~ — > 77 which is followed by the fragmentation of photons into vector mesons Vi(pi), V2(j>2)- Below these diagrams 
will be referred to as fragmentation diagrams. The second part are the diagrams shown in Fig. [TJ:, d will be referred 
to as non-fragmentation diagrams. The cross section cr(s) of the process e + e~ — > Vi(pi)p2(P2) can be written as the 
sum 

a ( s ) = + O-mt(s) + <Jnfr{s), (1) 

where cr/ r (s) and <y n t r {s) are the contributions due to the fragmentation and non-fragmentation diagrams correspond- 
ingly, <Ji n t{s) is the contribution of the interference between the fragmentation and non-fragmentation diagrams. 

In this paper double vector charmonium meson production (Vi = J/^,t/j') at B-factories will be considered. Com- 
monly, to study charmonium production one uses NRQCD formalism [lOj. In the framework of NRQCD charmonium 
mesons are considered as nonrelativistic systems with characteristic velocity v 2 ~ 0.3 and the amplitude of charmo- 
nium production is the series in the small parameter v 2 . The study of double vector charmonium meson production 



in e+e~ annihilation within NRQCD was carried out in papers [24 l25l l26l 1271. l28l. l29l|. 

In this paper light cone formalism [12J will be applied to the study of double vector charmonium meson production. 
Within this formalism the cross section is the series over inverse powers of characteristic energy of the process. In 
particular, at the energy </s — 10.6 GeV the expansion parameter for the process under consideration is 4M 2 /s ~ 0.4. 

To begin with let us determine the asymptotic behaviors of af r (s), tXj nt (s), a n f r (s) in the limit s — > oo. This can be 
done using the results of paper [281 ] - In Table [TJ the asymptotic behavior of <j/ r (s), <Ji nt (s), cr„/ r (s) in the limit s — > oo 
for different polarizations of vector mesons (Ai, A2) are shown. From this table one sees that at the leading order 
approximation in 1/s expansion only the fragmentation diagrams with the polarizations Xi = ±1, A2 = Tl contribute. 
It causes no difficulties to find the expression for this contribution 



<k±i,=Fi _ 167T 3 aV/i7f \/A / 1 - ri 



dx sM\M\ \{\ -a; 2 )A + 4rir 2 



(I-. 4 ), (2) 



where Mi,M^ are the masses of vector mesons, q c is the charge of c quark, x = cos 9, 9 is the angle between the 
momentums of electron and charmonium meson V\ , 

M 2 M 2 

ri = — ^, r 2 = —, A = l + rf + r|-2n-2r 2 -2rir 2 . (3) 
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Vi(Ai,pi)V r 2 (A 2 ,p2) 


a fr (s) 




Cn/r(s) 


Ai = ±1 A 2 = Tl 


~ 1/s 


~ 1/s 2 


~ 1/s 3 


Ai = ±1 A 2 = 


~l/s 2 


~ 1/s 3 


~l/s 4 


Ai = A 2 = ±1 








Ai = ±1 A 2 = ±1 


~l/s 3 


~l/s 4 


~l/s 5 


Ai = A 2 = 


~l/s 3 


~l/s 3 


~l/s 3 



TABLE I: The leading behavior of Uf r (s),a int (s),a n f r (s) in the limit s — > oo for the polarization of vector mesons Ai, A 2 . 



The constants f\ and /2 are defined through the matrix element of electromagnetic current J!f™ 

(y i ( ft ,A < )|JnO)=g c / i£ *(A i ). (4) 
This constants can be determined from the electronic width of vector meson Vi 

Formula ((2|) is valid for the production of different mesons. If two identical mesons are produced, this formula must be 
divided by 2. It should be noted here that formula (|2|) is in agreement with the result derived in paper [28|. To get the 
cross section of the process under consideration one should sum over all possible polarizations that give contribution 
to the cross section at the leading order approximation. Thus up to the corrections 0(1/ s 2 ) the cross section is 

ax ax \ s / 

In this section the cross section of the process e + e~ — ► Vi(pi)V^(p2) has been considered at the leading order 
approximation in 1/s expansion. Strictly speaking, to get the cross section at the leading order approximation one 
must expand formula in 1/s and keep only the first term. However, it turns out that the contribution of the 
fragmentation diagrams to the amplitude and the cross section can be calculated exactly. So, to improve the accuracy 
of the calculation done in here, the exact expression for cross section due to the fragmentation diagrams @ will be 
used. 

In the numerical calculation the following values of the parameters will be used: Mj/^ = 3.097 GeV, T(J/^ — > 
e+e~) = 5.55 ±0.14 KeV,Afy = 3.686 GeV,r(J/* -> e+e~) = 2.48 ±0.06 KeV 0. The results of the calculation 
are shown in Table ITT1 In the second colomn one can see the cross sections at the leading order approximation of 1/s 
expansion. In the third colomn the differential cross sections at the leading order approximation are integrated over 
the region | cos6*| < 0.8. 

There are different sources of uncertainty to the values of the cross section at the leading order approximation of 1/s 
expansion. The first one is QCD radiative corrections. These corrections can be divided into three groups. The first 
group is the radiative corrections due to the exchange of gluons between the quark and antiquark of one charmonium 
meson. Evidently, the same corrections appear as the radiative corrections to the electronic width of charmonium 
meson. These corrections are included into the values of the constants /j, which can determined from the electronic 
decay width of charmonia (O. So, the uncertainty due to the first group of the radiative corrections is reduced to the 
experimental uncertainty in the electronic decay widths of charmonia, which is few percents for J/ty and ip 1 mesons. 
The next group of the radiative corrections are the corrections due to the exchange of hard gluons between quarks 
or antiquarks of different mesons. This type of the corrections can be estimated as a 2 s (E)m 2 c / E 2 ~ 0.4% (28|, where 
E = y/s/2, m c is the mass of c quark. The last group of the radiative corrections is the corrections due to the exchange 
of soft gluons between quarks or antiquarks of different mesons, which can be estimated as (m c vfjE^ ~ 0.04% It 
is seen that the uncertainty due to the radiative corrections of the second and third group is very small. So, the main 
source of uncertainty can be reduced to the experimental uncertainty in the electronic decay widths of charmonium. 

Within light cone formalism in addition to QCD radiative corrections there are power corrections to the leading 
order approximation of 1/s expansion. These corrections appear due to the contribution of the fragmentation and 
non-fragmentation diagrams. In Table HT1 the error due to this corrections are estimated as M 2 /E 2 ~ 40%. To reduce 
this uncertainty let us consider 0(l/s 2 ) corrections to the leading order result. 
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III. NEXT-TO-LEADING ORDER CONTRIBUTION IN 1/s EXPANSION. 

To calculate the cross section at 0(1/ s 2 ) approximation of light cone formalism let us look to Table [U It is seen 
from this table that there are two contributions at this level of accuracy. The first one is due to the fragmentation 
diagrams with the following polarizations of the mesons A x = ±1, A 2 = and A x = 0, A 2 = ±1. It causes no 
difficulties to calculate these cross sections 

-±T = TMfMl ((l-x»)A + 4r 1 ^ ( (ri - ra + 1)a; + - ^ 2 + 4r 2 - 2) X + (n + r 2 - 1) ), 

df7 0,±i 16n 3 a 4 qjf 2 f 2 y/\ 2ri , 2 4 2 2 2 2 

-*T = ((l-^)A + 4^ 2 )^ (r2 " ri + 1)a: +(6-2-2^+4^-2), + [r x + r 2 - 1) ). (7) 

The second contribution arises from the interference between the fragmentation and non- fragmentation diagrams with 
polarization of vector mesons Ai = ±1, A 2 = =f1. It is not difficult to find the amplitude of the non- fragmentation 
diagrams (Fig. lc,d) for these polarization of vector mesons using LCF 

, f 2 9 7T 2 q 2 g gff 1 (^)ff 2 (^) r 

M = ^ L(2(p 2 fci) ~ 2(pifci))(e 1 e 2 )u(A;2>iu(A;i) + 

+s(e* 1 k 1 )u(k 2 )e^u(ki) + s(e 2 fci)M(fc 2 )eXfci))] A(cos 0, /z), (8) 

where w(fc 2 ),u(fci) are positron and electron bispinors, e^,e 2 are polarization vectors of charmonia. The constants 
gi(fj.) are defined as follows: 

(Viipi, X l )\Ca^C\0), l = gi(ji)(fr v - e* uP ^j. (9) 

It should be noted that the operator Ca a pC is not renormalization group invariant. For this reason the constant gi 
depends on scale as 

"M-dS)** *'- (10) 



The function A(x, (i) is defined as 

1 f j<- 0i(£ijM)<M6,M) 



A(x, l i) = ^J 



(i + 66) 2 -(a+6) 2 ^ 2 



xiyi X2V2 



(11) 



here Xi,a; 2 are the fractions of momentum carried by quark and antiquark in the first meson, yi,y 2 are the fractions 
of momentum carried by quark and antiquark in the second meson, £i = x\ — x 2l £ 2 = y\ — y 2 , </>i(£i, /i), ^2(^2, A*) are 
leading twist light cone distribution amplitudes of vector charmonium mesons with transverse polarization. 
Now some comments on formula © are in order: 

1. Formula © is the leading twist contribution to the amplitude of the diagrams shown in Fig[TJ;,d. For this reason 
it contains only the distribution amplitudes </>i(£i), ^2(^2) of the leading twist. 

2. It is seen from (J5J that the amplitude of the non- fragmentation diagrams depends on the distribution amplitudes 
4>i(£i, At) of vector mesons. If infinitely narrow distribution amplitudes & n) = S(£i) are substituted to formula 1(5)), 
than NRQCD result for the amplitude will be reproduced. If real distribution amplitudes </>i(£i,/z) are taken at scale 
/1 ~ m c , than formula ([8]) will resum the relativistic corrections to the cross section up to 0(1/ s 3 ) terms. To resum 
the relativistic and leading logarithmic corrections simultaneously one must take the distribution amplitudes <pi (£j , /j,) 
and the constants gi (/1) at the characteristic scale of the process /1 ~ ^fs. The calculation of the cross sections will be 
done at scale fi = E = y/s/2. 

The calculation of a±\ Tl will be done as follows. For the non-fragmentation diagrams the amplitude will be taken 
in form |8"]). For the fragmentation diagrams the exact expression for the amplitudes will be taken (see discussion in 
the previous section). Then the standard procedure for the calculation of the a±( Tl will be applied. Thus one gets 
the result 

dafl Tl _ 2 s ir 3 a i q*f 1 f 2 g 1 (E)g 2 (E) V\ ( 1 - n - r 2 
dx 3s 2 MfMf 

The total cross section has the form 



da da ! L x da± x dal r ±1 da+f^i ( 1 , 

— = 2 7' Tl +2— EM +2—^ +2 ELH + o( 13 

dx dx dx dx dx \s s 1 
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FIG. 2: The plots of the functions A(x, fi) for the processes e + e" -» J/* + J/* (fig. a), J/* + V' (fig. b), V' + V 1 ' (fig. c). 
Solid lines correspond to the leading order NRQCD predictions for the functions A(x,fi). Small dashed lines represent the 
functions A(x,fj,) if the relativistic corrections are taken into the account. Long dashed lines represent the functions A{x,fi) if 
the relativistic corrections and leading logarithmic radiative corrections are taken into the account simultaneously. 



IV. NUMERICAL RESULTS AND DISCUSSION. 



From formulas fT2")) , (flU)) one sees that the cross section depends on the function A(x,ii). This function takes 
into account internal motion of quark-antiquark pairs in mesons inside the hard part of the amplitude. In addition, 
this function resums the leading logarithmic radiative corrections to the amplitude. If one ignores both of these 
effects the function A(x,fi) equals unity, what corresponds to the leading order approximation of NRQCD. It is 
interesting to study how the relativistic and leading logarithmic radiative corrections can change the leading order 
NRQCD predictions. To do this one needs to know the distribution amplitudes <&(3lu) of IS and 2S states vector 
charmonium mesons. These distribution amplitudes were studied in papers [18L 1 19L 120J. The calculation of the 
functions A(x, fx) and the cross sections of the processes considered will be done using the models of distribution 
amplitudes proposed in these papers: 



M£,M~M C ) ~ (1-C 2 ) Exp 



•'25 



i-e 



(£,M~M C ) ~ (l-e)(a + e) Exp 



i-e 



(14) 



where M c = 1.2 GeV is the QCD sum rules mass parameter. For IS* charmonium state the constant j3 can vary within 
the interval 3.8 ± 0.7. For 2S charmonium state the constants a and (3 can vary within the intervals 0.031q'q 2 and 
2.5^0 8 correspondingly. 

Having models of distribution amplitudes (jT4j) , it causes no difficulties to calculate the functions A(x, fj,). The plots 
of the functions A(x, /i) for the processes e + e~ — > J + J/ty, J + ip' , ip' + ijj' are shown in Fig. [2^, b, c. Solid 
lines correspond to the leading order NRQCD predictions for the functions A(x,ii). Small dashed lines represent 
the functions A(x,[i) if relativistic corrections are taken into the account. Long dashed lines represent the functions 
A(x, /j,) if the relativistic corrections and leading logarithmic radiative corrections are are taken into the account 
simultaneously. 

From Fig. [2] one sees that the relativistic and leading logarithmic radiative corrections not only can change charac- 
teristic value of the function A(x, /i) but they can also considerably modify the shape of this function. This statement 
is especially true for the production of higher charmonia such as tp' meson, since the relativistic corrections play very 
important role in this case. 

Now let us calculate the cross sections of the processes under consideration. To do this one needs the values of 
the constants gi(fJ-) at some scale. Unfortunately, it is rather difficult to determine these constants directly from the 
experiment. The values of these constants can be obtained in the framework of NRQCD (see Appendix): 



gi(Mj /iS! ) =0.144 ±0.016 GeV 2 , g\ {M J/9 ) = 0.068 ± 0.022 GeV 2 



(15) 



The plots of the differential cross sections da/dx for the processes e + e~ — * J/^f + J/^, J/vP + i/j', ip' + ip' are shown 
in Fig. [3k, b, c. The values of the cross sections for the processes under study are shown in Table HT1 

There are different sources of uncertainty to the results obtained in this paper. The uncertainties in the fragmenta- 
tion contribution were discussed above. The uncertainties in the non-fragmentation contribution can be divided into 
the following groups: 
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do/tlx (fb) da/dx (fb) da/dx (fb) 




FIG. 3: The plots of the differential cross sections da/dx (x = cos 8) for the processes e + e~ — > J/* + J/\P (fig. a), J/ty + 
^' (fig. b), ip' + ip' (fig. c). Solid lines correspond to the 0(1/ s) contributions to the cross sections. Small dashed lines 
represent the cross sections at 0(1/ s 2 ) approximation. 



1. The uncertainty in the models of the distribution amplitudes 4>i(x,fi), which can be estimated through the 
variation of the parameters of these models (114|) . Thus it is not difficult to show that the error in the cross sections 
due to the uncertainty in the model of the distribution amplitude of J/^P meson is not very important (about few 
percents) and it will be ignored further. The error due the uncertainty in the model of the distribution amplitude 
of ip' meson is about 10% of the interference contribution for the process e + e~ — > J /\P + ip' and about 20% of the 
interference contribution for the process e + e~ — > ip' + ip' . 

2. The uncertainty due to the radiative corrections to the non-fragmentation diagrams. In the approach applied 
in this paper the leading logarithmic radiative corrections to the amplitude have been resummed in the distribution 
amplitudes. This fact allows us to estimate the rest of the radiative corrections as ~ a s (E) ~ 20%. It should be 
noted that if one does not resum the leading logarithmic radiative corrections the error of the calculation must be 
estimated as ~ a s (E) log s/Mj,^ ~ 50% instead of ~ a s (E) ~ 20% as it was done in paper [28[. 

3. The uncertainty due to the power corrections. This uncertainty is determined by the 0(1/ s 3 ) terms. One can 
estimate this source of uncertainty as ~ M 2 /E 2 ~ 40% 



0. 



4. The uncertainty in the values of constants (fTBj) . 

Adding all these uncertainties in quadrature one gets the total error of the calculations. 

Now it is interesting to compare the results for the cross sections with experimental data. The cross sections of the 
processes considered in this paper were measured at Belle collaboration [2j . Unfortunately, only the upper bound on 
these cross sections were determined: 

a(e + e- -► J/*J/*) x Br >2 (J/V) < 9.1 fb 90 % CL, 

a(e + e- -► J/Vip') x Br >2 (ip') < 13.3 fb 90 % CL, 

a(e + e- -> J/Vip') x Br >a (J/V) < 16.9 fb 90 % CL, 

a(e + e~ -► ip'ip') x Br >0 (ip') < 5.2 fb 90 % CL, (16) 

where -Br>2(U) denotes the branching fraction of V into final states with more than two charged tracks, Br>o(V r ) is 
the branching fraction of V into final states containing charged tracks. Unfortunately, the values of the -Br>o,2(U) 
are unknown. However, one can expect that the values of the Br>o(J/^), Br>2(J/4 r ), Bry^ip') are rather close 
to unity, what allows us to estimate a(e + e~ — * J/'fJ/'f) < 9.1 fb, a(e + e~ — > J/^ip 1 ) < 16.9 fb and a(e + e~ — > 
ip'ip') x Sr>o('0 / ) < 5.2 fb. These estimations are in agreement with the values of the cross sections obtained in this 
paper. 

From the results shown in Fig. [3] and in Table ILT1 one sees that 0(1/ s 1 ) contribution does not change greatly LO 
results for the process e + e~ — > J/^fJ/^f. The smallness of 0(1/ s 2 ) contribution for this process can be explained as 
follows. At 0(1/ s 2 ) approximation there are two contributions to the cross sections: the fragmentation diagrams and 
the interference of the fragmentation and non- fragmentation diagrams. These contributions are very near to each other 
and have different signs. So, due to partial cancellation 0(1/ s 2 ) contribution to the cross section is suppressed. If one 
further considers the production of higher charmonia, the value of the wave functions at the origin for these states 
are smaller and the total cross section becomes smaller. However, due to the relativistic and radiative corrections, 
collected in the functions A, the contribution of the non- fragmentation diagrams is enhanced (see Fig. [2]). For this 
reason 0(1/ s 2 ) contribution plays more significant role for the processes e + e~ — > J/^ip' and it is very important for 
the process e + e~ — * ip'ip'. 

It should be noted here that light cone formalism can be applied to study the production of light meson, for instance, 
p mesons. In this case, all formulas derived in this paper remain valid. The 0(1/ s 2 ) contribution can be estimated as 
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Vi v 2 


a LO (fb) 


cr |^?se|<o.8(f 3 ) 


a NLO (fb) 


„NLO m\ 

a \coae\<o.a\ ID ) 


a!2S](fb) 


l28 l rrw\ 


^(fb) 


J/# J/* 


2.12 ±0.85 


1.02 ±0.41 


2.02 ±0.25 


0.86 ±0.17 


1.69 ±0.35 


0.60 ±0.24 


1.8 - 2.3 




1.43 ±0.57 


0.77 ±0.31 


1.32 ± 0.16 


0.61 ±0.16 


0.95 ±0.36 


0.33 ±0.24 




^' v' 


0.24 ±0.10 


0.14 ±0.06 


0.20 ± 0.06 


0.10 ±0.05 


0.11 ±0.09 


0.04 ± 0.06 





TABLE II: The cross sections of the processes e + e~ — > J/tyJ/fy, J/^tj}', ip'ifj' . The second column contains the cross sections 
at the leading order approximation of 1/s expansion. The third column contains the differential cross sections integrated over 
the region | cos#| < 0.8. The values of the cross sections at 0(1/ s 2 ) approximation are shown in the forth and fifth columns. 
The sixth and seventh columns contain the results obtained in paper 28j. The results obtained in paper [29( ] are shown in the 
last column. 

~ M 2 /s, which is very small value for light mesons. For this reason, one can state that the values of the cross sections 
of double li ght meson production obtained in the approximation when only fragmentation diagrams are taken into 
account (27l.l28j are rather reliable. 

At the end of this section it is interesting to compare the results obtained in this paper with the results obtained in 
other papers devoted to the calculation of the same processes. It has already been noted that the processes of double 
vector mesons production were considered in the following papers 0, IH, [26|, [2?], [H, [29[ . In papers [H, [27| only the 
contribution arising from the fragmentation diagrams was considered. The results of their calculation are in good 
agreement with results obtained at LO approximation. The papers [24], |25l |28| , were written by the same group of 
authors, so it is reasonable to consider the results obtained in the last one [28[. In this paper the calculation was done 
at 0(v 2 ) approximation of NRQCD. The results obtained in paper [28| are shown in the last two columns of Table 
ITT1 It is seen that within the error of the calculation this results are in agreement with that obtained in this paper. 

Now let us consider the results obtained in paper [29j. In this paper the radiative corrections to the process 
e + e~ — > J/ty J/vE' were calculated. The result of the calculation can be written in the following form 

a 1 = a°(l ± ^K), (17) 

where a 1 is the cross section with the account of radiative corrections, a is the cross section without radiative 
corrections, the factor K = —11.19 for the pole mass of c-quark equals to 1.5 GeV. From this one sees that the 
radiative corrections to the cross section are very large, what leads to sizable reduction of the cross section a . As 
concerns the results obtained in this paper, it is seen from Tab. |TT] and from Fig. [21 [3] that the leading logarithmic 
radiative corrections do not change the cross section greatly. To the first sight, one can think that this contradicts 
to results [29| . However, there is no contradiction between these results. To see this let us consider the results of 
paper [29] in more detail. One of the input parameter for the calculation of the er° is the wave function at the origin 

|i? s J/ *(0)| 2 . In paper [291 ] this parameter was determined from the electron decay width of J/^f through the following 
formula 




It seen that this formula determines the value of \R S (0)| 2 taking into the account a s correction. As the result some 
part of the radiative corrections is present in a , which according to definition is the leading order in a s quantity. 
So, the authors of this paper separated the whole radiative corrections into two parts which nearly coincide but have 
different sign. If one now merges these two parts, the following result can be obtained 

a 1 = <t°(1 ± + (-11-19)—) = <t°(1 + (-0-52)^), (19) 

3 IT 7T 7T 

what is in agreement with the results obtained in this paper. The value of the cross section of the process e + e~ — ► 
J/tyJ/ty obtained in paper (2f| is presented in Table UTl The variation of the cross section is due to the variation of 
the pole mass of c quark 1.4 — 1.5 GeV. 
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APPENDIX A: THE CALCULATION OF THE CONSTANTS 

To calculate the values of the constants gi(/J.) ([9]) one can apply NRQCD formalism. At 0(v 2 ) approximation of 
NRQCD the constants gi(n) and can be written as follows jlO. [31j 

ff = (F i (e)| X +(ae)^|0)<0|x + (ae) ¥? |F i ( e ))x (i - ^ - l( v ^ , 

M = <^(£)|x + (a6>|0)<0| X + (a6>|^( e )) x ( 1 -^-l^l 0g ^--l{v 2 )X (Al) 



3 TT 3 TT 



where 



2 1 (0|x + (^(D)V|^(e)) 

^ (Olx+^^IViW) ' 1 j 

The calculation of the constants <7i(/x) will be done at scale /i = Mj/^. To diminish the error of the calculation let us 
consider the ratio gf(Mj/^)/ ff. At the same level of accuracy it can be written as follows 

The values of the constants gj(M//^) will be calculated with the following set of parameters: a s (Mj/ 9 ) = 0.25, 
( l;2 ),//* = 0-25 [12], (v 2 )^ = 0.54 [2(|. To estimate the error of the calculation one should take into account that 
within NRQCD the constant is double series in relativistic and radiative corrections. At NNLO approximation one has 
relativistic corrections ~ (v 2 ) 2 , radiative corrections to the short distance coefficient of the operator (0\x + (<7e)ip\Vi(e)) 

~ a 2 and radiative corrections to the short distance coefficient of the operator (0|x + (<?e)(D) 2 </?|14(e)) that can be 
estimated as ~ a s (v 2 ). Adding all these uncertainties in quadrature one can estimate the error of the calculation. 
Thus one gets 

gl(Mj^) = 0.144 ± 0.016 GeV 2 , 

al{Mj/*) = 0.068 ± 0.022 GeV 2 . (A4) 
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